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Abstract 

Die Verteilung von WasserstofFgas sowie die Rotationskurve der Milchstraße deuten auf die Existenz 
eines Rings aus kalter Dunkler Materie innerhalb der Milchstraßenebene in einem Abstand von rund 
14 Kiloparsec um das galaktische Zentrum. Es wird mit Hilfe eines semianalytischen Modells die 
Entstehung eines solchen Rings durch den Einfall von Zwerggalaxien in die Milchstraße simuliert. 
Die Entstehung des Rings kann sowohl durch einen Einzeleinfall als auch durch den Einfall mehrerer 
Zwerggalaxien erklärt werden, solange die numerische Exzentrizität der Orbits nicht 0.2 übersteigt 
und die Inklination der Orbits klein ist. Auch ein die Mindestmasse einer an der Ringentstehung 
beteiligten Zwerggalaxie beschreibender Zusammenhang ist ermittelt worden. 
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1 Einleitung 

Die unter Betrachtung der direkt beobachtbaren Materie 
zu hohen Geschwindigkeiten von Sternen in Galaxien und 
von Galaxien in Galaxienhaufen führten zur Annahme, 
dass sich in Galaxien und Galaxienhaufen mehr als nur 
die direkt beobachtbare Materie befi nden muss ( Zwickvl 



1937t lBosmalll98ll : Rubin et al.l . ll985 ). Diese Materie wird 



als Dunkle Material] bezeichnet, die Komponente mit Ge 
schwindigkeiten w ^ c als kalte Dunkle Materie (CDM). 
Der Umstand, dass CDM nicht direkt beobachtet werden 
kann, wird dadurch erklärt, dass sie nicht elektromagne- 
tisch wechselwirkt und dementsprechend keine elektroma- 
gnetischen Wellen aussendet. Diese Eigenschaft prägte den 
Namen der nur schwach wechselwirkenden, aber schweren 
Wimps (Weakly interacting massive particles). Das Stan- 
dardmodell kennt keine Teilchen, die diesen Anforderun- 
gen entsprechen, es existiert jedoch eine Vielzahl an Theo- 
rien, welche die Existenz solcher Teilchen postulieren. Die 
wohl populärste dieser Theorien ist die der Supersymme- 
trie. 

Aus der Verteilung von WasserstofFgas in der Milch- 
straße wurde auf die Massenverteilung in der Milchstra- 
ße geschlossen. In einem Abstand von 14kpc zum galak- 
tischen Zentrum liegt diesen Daten nach eine ringförmi- 
ge Ansammlung kalter Dunkler Materie mit einer Mas - 
se von bis zu 2.8 • 10^° Mq vor (jKalberla et al.l . l2007t ). 



Auch die Rotationskurve der Milchstraße deutet auf eine 
Massenansammlung bei R = 14kpc, siehe zum Beispiel 
Honma fc Sofud (1997a). Zusätzlich wird bei R = 14kpc 
eine Ansammlung von alten Sternen mit sich von an- 
deren Milchstraßensternen unterscheidender Geschwindig- 
keitsverteilung beobachtet (jlbata et al.l . 120031 ). Eine häufig 
diskutierte mögliche Ursache für die Entstehung dieses 
CDM-Rings bei R = 14kpc der Masse 2.8 • 10^° Mq ist, 
entsprechend dem Modell der hierarchischen Galaxienent- 
steh ung, der Einfall ein er einzelnen, massiven Zwerggala- 



xie (jCrane et al.l . l2003l ) 



Das Ziel dieser Arbeit ist die Entwicklung eines semi- 
analytischen Computermodells, welches die Simulation der 
Entstehung des äußeren CDM-Rings durch den Einfall von 
Zwerggalaxien ermöglicht und nutzt, um Aussagen über 
physikalische Eigenschaften der Zwerggalaxien und deren 
Orbits machen zu können. 

1.1 Modellbeschreibung 

Semianalytische Modelle haben n- Teilchen-Simulationen 
gegenüber den Vorteil, dass das Verständnis der funktiona- 
len Zusammenhänge des simulierten Systems leichter fällt. 
Um realistische und aussagekräftige Ergebnisse zu erhal- 
ten, müssen alle wesentlichen Prozesse des Modells einzeln 
beschrieben und in ihrer Wirkung verknüpft werden. Im 
Folgenden wird das dieser Arbeit zu Grunde liegende Mo- 
dell erläutert. 

Eine Zwerggalaxie befindet sich bei einem Radius zur 



deutlich schwereren Milchstraße. Die Zwerggalaxie bewegt 
sich durch ein Gebiet der Dichte p im Gravitationspotenti- 
al der Milchstraße mit der Anfangsgeschwindigkeit Vi. Die 
Modellierung der Dichteverteilung der Milchstraße ist in 
Abschnitt 12.11 die des Gravitationspotentials in Abschnitt 
12.21 beschrieben. Bei der Bewegung im Gravitationspoten- 
tial der Milchstraße verliert die Zwerggalaxie durch dyna- 
mische Reibung Drehimpuls und Bewegungsenergie. Diese 
Effekte der dynamischen Reibung im Halo werden in Ab- 
schnitt 12.31 behandelt. Bei jedem Durchgang der Zwerg- 
galaxie durch die galaktische Scheibe kommt es zu Kolli- 
sionen zwischen den HI-Atomen des Wasserstoffgases der 
Zwerggalaxie und dem Wasserstoffgas innerhalb der galak- 
tischen Scheibe der Milchstraße. Dadurch wird das Was- 
serstoffgas der Zwerggalaxie abgebremst und sammelt sich 
in der galaktischen Ebene der Milchstraße. Die Kollisio- 
nen werden in Abschnitt 12 .41 beschrieben. Durch die dyna- 
mische Reibung spiralisiert die Zwerggalaxie in Richtung 
des galaktischen Zentrums der Milchstraße. Dabei steigen 
mit sinkendem Abstand zum galaktischen Zentrum die auf 
die Zwerggalaxie wirkenden Gezeitenkräfte, die bei einem 
Abstand rt zum galaktischen Zentrum zum Zerreißen der 
Zwerggalaxie führen. Auf die Gezeitenkräfte wird in Ab- 
schnitt 12.51 näher eingegangen. Die zerrissene Zwerggala- 
xie bildet einen Sternstrom (sowohl aus Sternen als auch 
aus CDM und Wasserstoffgas bestehend) , der sich zu einer 
ringförmigen Struktur um das galaktische Zentrum entwi- 
ckelt. Dieser Sternstrom ist aufgrund der ringförmig und 
weitläufig verteilten Masse unanfällig für dynamische Rei- 
bung des Kalos, ändert daher seinen Abstand zum galak- 
tischen Zentrum nur noch wenig. Bei dem Einfall einer 
Zwerggalaxie mit geringer Inklination kann sich also im 
Laufe der Zeit eine innerhalb der galaktischen Ebene lie- 
gende, ringförmige Struktur aus Sternen und CDM bilden. 

2 Simulation 

In diesem Abschnitt wird die Modellierung der einzelnen 
in Abschnitt 11.11 erwähnten Prozesse beschrieben. Des wei- 
teren wird das verwendete Verfahren zur Zeitintegration 
der Flugbahn in Abschnitt 12.61 vorgestellt. Die Implemen- 
tierung der Simulation erfolgt objektorientiert in der Pro- 
grammiersprache cH — h, auf diese wird allerdings in dieser 
Arbeit nicht näher eingegangen. 

2.1 Dichteverteilung in der Milchstraße 

Das verwendete Modell der Milchstraße lässt sich in drei 
Komponenten aufteilen: Die galaktische Scheibe und den 
zentralen Bulge, welche überwiegend aus Wasserstoffgas 
bestehen, sowie das Halo aus Dunkler Materie. 

Die Dichteverteilung (r) des sphärisch sym metrischen 
Halos wird mit Hilfe des von IS aha et al. Iläooi) behandel- 
ten Dichtepr ofils beschrieben, we lches für r ^ dem 
NFW-Profil (jNavarro et al.l . Il996l) entspricht: 



^Ein Überblick über die Erforschung der Dunklen Materie bis 1999 
ist von S. van den Bergh gegeben, zu finden bei arXiv:astro- 
ph/9904251vl 



Phir) = Pom ■ 



1 



(1) 
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Die ebenfalls sphärische Dichtever teilung Pb{r) des Bulges 
wird wie von iDehnen et al.l (|l998l ). allerdings ohne Exzen- 
trizitäten, modelliert: 



Pb{r) = Pom- (—) exp(^) 



(2) 



Für r ^ rt^fc ist damit pb{r) proportional zu r '''' und 
nimmt sobald r > r^^ stark ab. Die Dichte po.fci der Radius 
ro.fc sowie der Exponent jb sind durch Beobachtungsdaten 
des COBE/DIRBE Satelliten bestimmt worden. 

Die Dichte der gala ktischen Schei b e (Dis k) wird eben- 
falls nach einem von Dehnen et al. (|l998l ) vorgestellten 
Modell beschrieben. In diesem Modell besteht die galak- 
tische Scheibe aus dem Interstellaren Medium (ISM) so- 
wie der dünnen und der dicken Sternscheibe. Die folgende 
Gleichung beschreibt die Dichte der Disk: 



Prf(r, z) = pa,d ■ exp 



R 



7n.d 



Rd 



(3) 



Dabei sind Rm,d und Rd Radien zur Skalierung der Schei- 
be, die aus Beobachtungsdaten g ewonnen wurden. Eben- 



so sind von iDehnen et al.l (jl998D die Ausdehnungen a 



bezüglich der galaktischen Ebene sowie die Dichten po.d 
für die drei Komponenten der Scheibe ermittelt worden. 

Die Gesamtdichte pg der Milchstraße berechnet sich aus 
der Summe der Dichten von Halo, Bulge und Disk: 



Pg{r, z) = ph{r) + pb{r) + pd{r, z) 



(4) 



Vereinfachend wird die Flächendichte PR{r) des durch 
den Einfall von Zwerggalaxien gebildeten Rings aus CDM 
der Masse über die Breite Rmax — Rmin des Rings in 
diesem Modell als normalverteilt angenommen, der mittle- 
re Radius des Rings R stellt als Erwartungswert der Dich- 
tefunktion den Radius höchster Dichte dar. Die Standard- 
abweichung an beschreibt die Größe des Bereichs hoher 
Dichte im CDM- Ring und berechnet sich über d = 2a r, 
wobei d den Durchmesser der verursachenden Zwerggala- 
xie darstellt. 



1 



Mr. 



2nR crßV27r 



exp 



{R- 



2 4 



(5) 



In Tabelle [T] sind die Parameter der verwendeten Model- 
le aufgelistet und Abbildung [T] zeigt die Dichteverteilung 
in der Ebene der Milchstraße {z = 0) in Abhängigkeit vom 
Abstand r zum galaktischen Zentrum. 

2.2 Gravitationspotential 

Die Orbits von in die Milchstraße einfallenden Zwergga- 
laxien hängen von ihrer Gesamtenergie Eges ab, welche 
sich aus der kinetischen Energie E^in und der potentiellen 
Energie Epot zusammen setzt. Die potentielle Energie Epot 
lässt sich aus dem Produkt der Masse der Zwerggalaxie 
m und dem Gravitationspotential $(r) berechnen. Dabei 
hängt <I'(r) vom Positionsvektor r der Zwerggalaxie sowie 
von der in Abschnitt 12.11 beschriebenen Dichteverteilung 
p{r) ab. Für sphärisch symmetrische Dichteverteilungen 
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Tabelle 1: Param eter der Dichteverteilung en gemäß 
iDehnen et äP l| 19981 ) (Disk und Bulge) sowie ISaha et al.l 
(|2009l 'l (Halö). Dichten in Mo/kpc^, Strecken in kpc 
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Abbildung 1: Dichteverteilung p(r) in der galatischen Ebene der 
Milchstraße in Abhängigkeit vom Abstand r zum galaktischen 
Zentrum 



wie die des Halos und die des Bulges ist ^sph nur von 
r = |r| abhängig: 



GM(s) 



ds mit M(s) ^ 4tt a p{a)da (6) 



Dabei beschreibt M(r) die Masse innerhalb der Kugel mit 
dem Radius r, welche bei sphärischen Massenverteilungen 
gemäß dem newtonschen Theorem als Punktmasse bei r = 
betrachtet werden kann. 

Für ringförmige Massenverteilungen innerhalb der ga- 
laktischen Ebene um das Zentrum der Milchstraße ist das 
Gravitationspotential in einem Punkt P abhängig von der 
Position des Punktes bezüglich der galaktischen Ebene, 
siehe Abb. [2] Der Abstand D/j vom Punkt P zu einem 




Abbildung 2: Ring um das galaktische Zentrum mit Radius R 
und Masse Mu innerhalb der galaktischen Ebene. Der Schnitt- 
punkt der gestrichelten Linien beschreibt die Lage des galakti- 
schen Zentrums. 

Punkt auf dem Ring mit Radius R, abhängig vom Winkel 
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6, berechnet sich wie folgt: 

Dfl(r, z, d) = (r^ + + ~ 2Rr cos{d) 



1/2 



(7) 



Der CDM-Ring wird als Ring der Breite Rmax — Rmin 
innerhalb der galaktischen Ebene modelliert. Es ergibt 
sich für das Gravitationspotential des Rings ^ring — 
—Gdm/ Du durch Integration über alle Radien R des 
Rings und den Winkel 5: 



DR{r,z,S) 



dS dR 



(8) 



Dabei beschreibt pj^ die Masse pro Umfang des Rings bei 
einem bestimmten Radius R. Um Rechenzeit zu sparen, 
wird aufgrund der im Vergleich zum CDM-Halo kleinen 
Masse der Disk die Masse der Disk als Näherung in ihrer 
gravitativen Auswirkung über das Potential des sphärisch 
symmetrischen Halos beschrieben. 

2.3 Dynamische Reibung im Halo 

Bewegt sich eine Zwerggalaxie mit der Geschwindigkeit 
V durch ein Gebiet der Milchstraße der Dichte p{r), so 
übt diese eine Kraft proportional zu ihrer Masse m auf 
umliegende Körper, wie Sterne oder Gaswolken, aus. Die- 
se Körper erfahren eine Beschleunigung in Richtung der 
Zwerggalaxie und führen damit im Bereich in Bewegungs- 
richtung hinter der Zwerggalaxie zu einem Dichteanstieg. 
Durch die Masse in diesem Bereich höherer Dichte wirkt 
nun eine Kraft auf die Zwerggalaxie, die entgegen ihrer 
Bewegungsrichtung gerichtet ist. Diese Kraft führt zu ei- 
ner Abbremsung der Zwerggalaxie, welche daher Bewe- 
gungsenergie und Drehimpuls verliert. Ers tmalig beschrie- 
ben w urde die dynamische Reibung von IChandrasekhar 
( 1943h . Der Prozess ist besonders effektiv bei massiven, 
langsamen Zwerggalaxien, da diese für eine lange Zeit eine 
große Kraft auf umliegende Körper ausüben, welche eine 
große Kraft entgegen der Bewegungsrichtung der Zwerg- 
galaxie bewirken. 

In den äußeren Bereichen des sphärischen Halos der 
Milchstraße können die Neigungen der Orbits von Teil- 
chen im Halo als zufällig verteilt angesehen werden. In die- 
sem Bereich hat sich, aufgrund der zufälligen Orbitneigun- 
gen, die Modellierung der Geschwindigkeiten d er Teilchen 
des H alos über eine Maxwellverteilung bewährt ( Tremain3 . 
1976h . Unter diesen Annahmen lässt sich mit der folgenden 
Gleichung die Änderung i) pro Zeit der Geschwindigkeit v 
der Zwerggalaxie bestimmen: 



-47rG^m(?(i?) InA— 



''erf( — ) — —eri' ( — 1 

Vr Vr. Vr 



(9) 



Dabei ist G die Gravitationskonstante und InA ist der 
Coulomb-Logarithmus, welcher die Reichweite der dyna- 
mischen Reibung festlegt. Dieser wird über eine Compu- 
tersimulation abgeschätzt: Ein Probekörper bewegt sich 
durch einen das Milchstraßenhalo repräsentierenden Raum 
mit 1000 Massepunkten, die gravitativ mit dem Pro- 
bekörper wechselwirken und ihn durch dynamische Rei- 



bung abbremsen. Aus der Simulation wird die Geschwin- 
digkeitsänderung pro Zeit v bestimmt und über Glei- 
chung |Ü] auf den entsprechenden Coulomb-Logarithmus 
InA geschlossen. Die Simulation führt bei einem Massen- 
feld der Dichte po./i und einem Probekörper der Masse 
m = 2.8 • 10^" M0 und der Geschwindigkeit — 250km/s 
zum Coulomb-Log arithmus InA = 3.01. Die ser Wert liegt 
im Bereich der von I Just fc Penarrubial (|2005[ ) verwendeten 
Coulomb-Logarithmen für das Milchstraßenhalo. 



5 
4 
3 
2 

1 ^ 



•S8 Ja- < 




Abbildung 3: Probekörper (ausgefüllte Kugel) im Massenfeld, 
wobei die Bewegungsrichtung des Probekörpers der Richtung 
der positiven y-Achse entspricht. In dem (lOkpc)'^ großen Vo- 
lumen befinden sich 1000 Massenpunkte. 

Die Simulationen sind auf numerische Konvergenz getes- 
tet. Dazu wurde in einem Simulationslauf die Schrittweite 
A< zur numerischen Integration der Bewegung der Partikel 
um Faktor zehn verringert. In einem anderen Simulations- 
lauf wurde die Anzahl der Massenpunkte um Faktor zehn 
erhöht. Der jeweils ermittelte Coulomb-Logarithmus InA 
weicht um weniger als 1 Prozent von InA = 3.01 ab. 

2.4 Wechselwirkung von Baryonen und Disk 

Während die Dunkle Materie beim Einfall der Zwerggala- 
xie nur über dynamische Reibung Bewegungsenergie ver- 
liert, kann die baryonische Materie aus der Zwerggalaxie 
durch Kollision mit der baryonischen Materie der Milch- 
straße abgebremst werden. Da Wasserstoff das Element 
mit der größten Häufigkeit in der Milchstraße ist und stark 
konzentriert in der Disk vorliegt, werden die Kollisionen als 
solche zwischen HI-Atomen simuliert. Die Stöße werden 
als elastisch und zwischen Kugeln des Durchmessers d mit 
harten Schalen beschrieben. Ziel ist es in diesem Abschnitt, 
eine Durchschnittsgeschwindigkeit eines HI-Atoms nach 
einer Kollision zu berechnen. 

Die thermische Bewegung der Hl-Atomc von Zwerggala- 
xie und Milchstraßenscheibe führt, durch die geringe Dich- 
te des Wasserstoffgases und damit durch die große freie 
Weglänge zwischen den HI-Atomen, zu Geschwindigkeits- 
unterschieden zwischen den einzelnen HI-Atomen. Die Ge- 
schwindigkeit eines Wasserstoffatoms in der folgenden Be- 
rechnung beschreibt stets die Durchschnittsgeschwindig- 
keit der Wasserstoffatome von Zwerggalaxie oder Milch- 



4 



Raphael Errani: Modell zur Entstehung eines Rings aus Dunkler Materie in der Milchstraßenebene 



straßenscheibe, die thermischen Geschwindigkeiten mit- 
tein sich damit zu null und werden in der Berechnung nicht 
beachtet. 

In der folgenden Betrachtung wird ein HI-Atom als ru- 
hend angenommen, das andere Atom bewegt sich mit der 
Relativgeschwindigkeit v. Die Richtung der x- Achse des 
genutzten Koordinatensystems entspricht jener von v, sie- 
he Abbildung Entscheidend für die Kollision ist der 



Stoßparametes b daher: 




Abbildung 4: Koordinatensystem, in welchem die Kollision be- 
trachtet wird. Die x-Achse ist so gewählt, dass ihre Richtung 
jener der Relativgeschwindigkeit v der HI-Atome entspricht. 

Kollisionwinkel 7. Bei einer Kollision verliert das stoßen- 
de Atom die Geschwindigkeitskomponente parallel zur 
Verbindungslinie der beiden HI-Atome an das ruhende 
Atom und behält selbst lediglich die Geschwindigkeitskom- 
ponente v± senkrecht zur Verbindungslinie. 

Der Winkel 7 tritt um die x- Achse des verwendeten Ko- 
ordinatensystems im Raum rotiert auf. Aus Symmetrie- 
gründen bleibt bei Betrachtung aller möglichen Rotatio- 
nen von 7 lediglich die x-Komponente v' von v± erhalten, 
da sich die y-Komponenten bei hinreichend großer Stoß- 
anzahl zu null mittein. Ist der Winkel 7 bekannt, so ergibt 
sich für v' (vergleiche Abbildung S]) : 



«'(7) = V sin^ 7 



(10) 



Im nächsten Schritt wird die Wahrscheinlichkeit eines 
Stoßwinkels 7 berechnet. Die Wahrscheinlichkeit einer 
Kollision hängt von der Größe des Stoßparameters b 
ab: Kleine Stoßparameter, die zu frontalen Kollisionen 
führen, haben kleine Wahrscheinlichkeiten; große Stoßpa- 
rameter, die lediglich das gegenseitige Streifen der HI- 
Atome bewirken, haben höhere Wahrscheinlichkeiten. Ab- 
bildung [Ü] zeigt ein Stoßparameterintervall von bi bis 
62. Die Wahrscheinlichkeit dieses Stoßparameterintervalls 
entspricht seinem Flächenverhältnis zum Wirkungsquer- 
schnitt a. Bei bekanntem Wirkungsquerschnitt a = ttcP 




A b 



o 
b2 



Abbildung 5: Wirkungsquerschnitt a und Stoßparameterinter- 
vall zwischen fei und fe2. Dem Radius d entspricht der Abstand 
der beiden HI-Atome zum Zeitpunkt der Kollision. 

gilt für die kumulierte Wahrscheinlichkeit Pcum{b) eines 



[TTb^lb 
- 



(11) 



Aus Abbildung |4] ist zu erkennen, dass b = dsiwy. Setzt 
man dies in Gleichung [TT] ein und leitet ab, um zur partiel- 
len Wahrscheinlichkeit eines Winkels 7 zu gelangen, erhält 



d sin^ 7 . 
Pparh) = — -j^ — = 2sm7COS7 



(12) 



Da nun die Wahrscheinlichkeit eines Stoßwinkels 7 sowie 
die mit einem bestimmten 7 resultierende Geschwindig- 
keit v' bekannt ist, lässt sich über Gewichtung von v'{'y) 
mit Ppar{l) und Integration über alle möglichen 7 auf 
die durchschnittliche Geschwindigkeit v' nach einem Stoß 
schließen: 



v' ~ 2v f sm^ 7 cos 7^7 — 



(13) 



Für die Geschwindigkeit Vtrans eines HI-Atoms des Was- 
serstoffgases der Zwerggalaxie nach Transit der Disk der 
Milchstraße muss die Flächendichte pni an HI-Atomen in 
der Disk bekannt sein. Aus l/pHi ergibt sich die durch- 
schnittliche jedes HI-Atom umgebende Fläche Ahi- Das 
Verhältnis von a und dieser Fläche Ahi beschreibt die 
Wahrscheinlichkeit für eine Kollision, siehe Abbildung [ü] 
Es ergibt sich aus diesen Überlegungen die Geschwindig- 



Q H1 




Abbildung 6: Fläche Ah\ um ein HI-Atom sowie Wirkungs- 
querschnitt (T. Die Wahrscheinlichkeit einer Kollision ist durch 
das Flächenverhältnis von a und Ahi gegeben. 

keit Vtrans cincs HI-Atoms nach Transit der galaktischen 
Scheibe über: 



I2J 



(14) 



Dabei ist v die Relativgeschwindigkeit des HI-Atoms der 
Zwerggalaxie zu den HI-Atomen der Disk bei Eintritt in 
die Disk. 

2.5 Gezeitenkräfte 

Fällt eine Zwerggalaxie in die Milchstraße ein, so steigen 
mit abnehmendem Abstand zum galaktischen Zentrum die 
auf die Zwerggalaxie wirkenden Gezeitenkräfte. Ist die Dif- 
ferenz Fm,i — Fm,2 = der auf die Zwerggalaxie wir- 
kenden Kräfte im Gravitationspotential der Milchstraße 
größer als jene die Zwerggalaxie zusammen haltende Kraft 
Fzw , wird die Zwerggalaxie zu einem Sternstrom zerrissen. 
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Die Kraft F^^ wird vereinfachend als jene Gravitations- 



M,1 



M,2 



Abbildung 7; Kräfte zwischen zwei Massen Mz„/2 im Gravita- 
tionspotential der Milchstraße 

kraft beschrieben, die zwischen zwei Körpern der Masse 
Mzwl2 im Abstand d wirkt, wobei M^w die Masse und 
d den Durchmesser der Zwerggalaxie darstellt. Um AF/\f 
zu berechnen, wird die Differenz der Kräfte zwischen der 
Milchstraße und der Masse M^^j^ in einem Abstand R 
und R + d zum galaktischen Zentrum bestimmt. Dabei ist 
R der Abstand des galaktischen Zentrums zur dem galak- 
tischen Zentrum zugewandten Seite der Zwerggalaxie. Es 
folgt nach einigem Umformen, dass die Zwerggalaxie stabil 
ist, so lange gilt: 



,^ ^ M{r) 



2d2 



M{r + d) 

{r + dy 



(15) 



Dabei ist M{r) die Masse der Milchstraße innerhalb des 
Radius r. Für Radien kleiner als rt ist die oben genannte 
Bedingung nicht mehr erfüllt und die Zwerggalaxie wird 
zu einem Sternstrom zerrissen. 

2.6 Numerische Integration 

Die Bewegung der Zwerggalaxie wird durch das Gravita- 
tionspotential $(r) der Milchstraße, siehe Abschnitt 12.21 
bestimmt. Die von der Zwerggalaxie erfahrene Beschleuni- 
gung f ergibt sich aus der Ableitung des Gravitationspo- 
tentials nach dem Positionsvekor und führt zur folgenden 
Differentialgleichung 2. Ordnung: 



r = — 



dr 



(16) 



Diese wird durch zweifaches Anwenden des vi erstufigen 
Runge Kutta- Verfahrens (siehe zum Beispiel iHellingsl . 
[l994) numerisch integriert. Um aus der Beschleunigung r 
auf die Geschwindigkeit r zu schließen, werden vier Stei- 
gungen kf benötigt: 



k^T 



(17) 
(18) 
(19) 
(20) 



Dabei ist At die für die Genauigkeit der Integration ent- 
scheidende Schrittweite. Die Geschwindigkeit r ergibt sich 
aus dem gewichteten Mittel der 4 Steigungen nach Multi- 
plikation mit der Schrittweite: 



■Tt+l 



kir + 2k2r + 2^3^ + k^j. 



At 



(21) 



Die Steigungen kf- müssen abwechselnd mit den Steigun- 
gen kr berechnet werden: 



kir = r 

K2-r = r + klr — 

. , At 

K'ir = r + k2r — 
kir ~ r + fcsrAt 

Der neue Positionsvektor r ergibt sich aus: 

kir + 2Ä:2r + Sfcsr + k^r 



rt+1 = rt 
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-At 



(22) 
(23) 

(24) 
(25) 

(26) 



Da r eine vektorielle Größe ist, muss die eben beschriebe- 
ne Zeitintegration für jede Komponente einzeln ausgeführt 
werden. 



3 Ergebnisse 

In diesem Abschnitt werden die Ergebnisse der in Ab- 
schnitt [5] beschriebenen Simulation dargestellt. 

Das sphärisch symmetrische Dichteprofil von Halo ph (r) 
und Bulge pb(r) wird mit einer räumlichen Auflösung von 
l/30kpc diskretisiert. Mit Gleichung |ü] und numerischer 
Berechnung von M{r) wird damit das sphärisch symme- 
trische Gravitationspotential beschrieben, aus Gleichung 
Owird auf das Gravitationspotential des Rings aus CDM 
geschlossen. 

Die Schrittweite At der Zeitintegration beträgt 1 • 10^ 
Jahre. Es wird ein Zeitraum dem Alter des Universums 



entsp rechend von 13.7-10^ Jahren simuliert (jSpergel et al 
20031) . 



3.1 Exzentrizität des Orbits 

Die bei 14 kpc Entfernung zum galaktischen Zentrum be- 
obachtete Ansammlung von kalter Dunkler Materie, Ster- 
nen und Wa sserstoffgas wird als a n nähernd rin g förmi g an- 



genommen (jKalberla et al.l . 120071 : llbata et al.l . l2003f ). so 



mit kann von einer kleinen numerischen Exzentrizität e 
des Rings ausgegangen werden. 

Es soll nun geklärt werden, in wie weit die Exzentrizität 
des Orbits der Zwerggalaxie sich auf den entstehenden 
Ring auswirkt. Die Exzentrizität des Orbits der Zwerg- 
galaxie lässt sich, da aufgrund der wirksamen dynami- 
schen Reibung und der Massenverteilung M(r) const die 
Zwerggalaxie keine reinen Keplerellipsen beschreibt, nicht 
über die herkömmliche numerische Exzentrizität einer El- 
lipse darstellen. In Anlehnung an die Exzentrizität einer 
Ellipse e wird hier die Exzentrizität St wie folgt definiert: 



Et - 



'^t,max H" '^t,min 



(27) 



Dabei entsprechen die Radien rt^max und rt^min dem je- 
weils zeitlich letzten maximalen und minimalen Abstand 
der Zwerggalaxie zum galaktischen Zentrum. 

Um das Verhalten von et in Abhängigkeit der Zeit t zu 
untersuchen, wird eine Zwerggalaxie beim Radius = 100 
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Abbildung 8: Orbits einer Zwerggalaxie mit der Masse m = 2.8- lO^'^M© im Gravitationspotential der Milchstraße. Es sind Orbits 
mit verschiedenen Exzentrizitäten d zu Beginn der Simulation dargestellt. Die Inklination ist bei allen Orbits mit i = gewählt. 
Die Galaxie startet beim Radius Vi = 100 kpc und bewegt sich in Richtung des galaktischen Zentrums (durchgängige Linie). 
Bei Erreichen von rt wird sie durch Gezeitenkräfte zerrissen. Das Orbit ist von dort an mit gestrichelter Linie weitergeführt, 
entlang d ieser Bahn werden sich Sterne, Wasserstoffgas und CDM der Zwerggalaxie ansammeln. Die Punkte markieren die von 
ICrane et al.. (,2003i 'l beobachteten Sterne der Ringstruktur. 



kpc mit der Geschwindigkeit Vi gestartet und bei jedem 
Erreichen des Perizenters die jeweilige Exzentrizität St be- 
stimmt. Dabei stellt den größten Abstand zum galak- 
tischen Zentrum dar. Für diesen Punkt wird für verschie- 
dene numerische Exzentrizitäten Ci einer Keplerellipse die 
jeweilige Ap oapsisgeschwindigkeit b erechnet (siehe dazu 
zum Beispiel iKarttunen et al. . 2007 ). welche als Anfangs- 
geschwindigkeit Vi = Vc\/1 — Ci der Zwerggalaxie verwen- 
det wird. Dabei ist Vc die stabile Kreisbahngeschwindigkeit 
für Vi. Abbildung [9] zeigt die Entwicklung der Exzentrizität 
e in Abhängigkeit der Zeit t für verschiedene numerische 
Exzentrizitäten e^, die zum Berechnen der Startgeschwin- 
digkeit verwendet wurden. 

Nach der Zerstörung der Zwerggalaxie durch Gezeiten- 
kräfte liegt diese als Sternstrom vor und erfährt aus diesem 
Grund weniger dynamische Reibung. Dieser Sternstrom 
behält daher die Exzentrizität des Orbits der ihn verur- 
sachenden Zwerggalaxie weitgehend bei und führt auch zu 
einem Ring entsprechend der Exzentrizität e des Orbits 
der Zwerggalaxie. 

Die Exzentrizität e des Orbits nimmt im Laufe der Zeit 
durch Effekte dynamischer Reibung ab, der Orbit wird al- 
so kreisförmiger. Das ist ein häufig beobachteter Prozess 



bei S imulationen mit konstantem Coulomb Logarithmus 
InA ijHashimotd . I20Q3I ). Dieser ist jedoch zu langsam, um 
Orbits mit numerischen Exzentrizitäten > 0.2 so weit 
kreisförmig zu machen, als dass die Masse der Zwerggala- 
xie in der beobachteten Weise um das galaktische Zentrum 
verteilt würde, siehe Abbildung [S] Daher kann die anfäng- 
liche numerische Exzentrizität des Orbits einer eingefal- 
lenen und an der Ringenstehung beteiligten Zwerggalaxie 
nicht größer als etwa 0.2 gewesen sein. 

3.2 Inklination des Orbits 

Kreuzt die Zwerggalaxie auf ihrem Orbit die galaktische 
Scheibe der Milchstraße, kommt es zu Kollisionen zwischen 
dem interstellaren Wasserstoffgas der Milchstraßenschei- 
be und dem Wasserstoffgas der Zwerggalaxie, siehe Ab- 
schnitt 12.41 Je dichter die Zwerggalaxie dabei dem galakti- 
schen Zentrum ist, desto höher ist die Dichte des Was- 
serstoffgases an der Stelle, an welcher die Zwerggalaxie 
die Milchstraßenscheibe durchkreuzt, siehe Gleichung [3] 
Damit steigt mit sinkendem Abstand zum galaktischen 
Zentrum die Wahrscheinlichkeit einer Kollision zwischen 
Wasserstoff aus der Milchstraße und Wasserstoff aus der 
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Abbildung 9: Entwicklung der Exzentrizität e des Orbits einer 
Zwerggalaxie in Abhängigkeit der Zeit t. Der Zusammenhang 
ist für unterschiedliche Werte der numerischen Exzentrizität 
am Anfang der Simulation dargestellt. Das Ende einer Linie 
entspricht dem Zeitpunkt der Zerstörung der Galaxie durch 
Gezeitenkräfte, nach welchem et nahezu konstant bleibt. 



Zwerggalaxie. Abbildung [TÜ] zeigt für drei Zerggalaxien der 
Masse 2.8 • IO^^Mq die Entwicklung der Inklination i im 
Laufe der Simulation für verschiedene anfängliche Inklina- 
tionen ii und verschiedene Startradien des Orbits, die 
anfängliche Exzentrizität beträgt 0.1. Die Abbildung 
zeigt, dass bereits nach wenigen Umläufen um das galakti- 
sche Zentrum die Bewegungsenergie der Geschwindigkeits- 
komponente der Zwerggalaxie senkrecht zur galaktischen 
Scheibe vom Wasserstoffgas der Milchstraße aufgenommen 
wurde. Dadurch sinkt die Inklination des WasserstofFga- 
ses der Zwerggalaxie schnell. Die vom Wasserstoffgas der 
Milchstraße aufgenommene Bewegungsenergie kann durch 
weitere Kollisionen innerhalb der galaktischen Scheibe in 
Wärmeenergie umgewandelt werden. 

Die Dunkle Materie wechselwirkt nicht elektromagne- 
tisch, siehe Abschnitt [T] Deshalb kann die Inklination der 
Dunklen Materie der Zwerggalaxie nicht wie die Inklinati- 
on des Wasserstoffgases durch Kollision abnehmen. Durch 
dynamische Reibung in abgeflachten CDM-Halos ist ein 
Prozess bekannt, welcher beschreibt, wie sich rein gravita- 
tiv beeinflusst di e Inklination des Orbits einer Zwerggala- 
xie ändern kann ( Penarrubia et al.l . 12004 ) . Dieser findet in 
dem hier verwendeten Modell jedoch keine Anwendung, da 
das simulierte CDM-Halo sphärisch symmetrisch ist, sie- 
he Gleichung [TJ Fällt eine Zwerggalaxie in die Milchstra- 
ße mit einem Orbit deutlich außerhalb der galaktischen 
Ebene ein, kommt es zu einer Trennung von Wasserstoff- 
gas und Dunkler Materie der Zwerggalaxie, da das Was- 
serstoffgas Bewegungsenergie an die Milchstraßenscheibe 
abgeben und dabei die Inklination verringern kann, die 
Dunkle Materie hingegen nicht. 



3.3 Masse der Galaxie 

Damit eine Zwerggalaxie vom Startradius r.i bis zum Radi- 
us rt , an welchem sie durch Gezeitenkräfte zerrissen wird, 
gelangt, muss sie über dynamische Reibung (siehe Ab- 
schnitt 12.3p Bewegungsenergie abgeben. Die in der simu- 
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Abbildung 10: Zeitliche Entwicklung des Abstands Vi der 
Zwerggalaxie zum galaktischen Zentrum (oberes Diagramm) 
sowie der Inklination i des interstellaren Wasserstoffgases der 
Zwerggalaxie beim mehrmaligen Durchkreuzen der galakti- 
schen Scheibe (unteres Diagramm). 



Herten Zeit über dynamische Reibung maximal abzuge- 
bene Bewegungsenergie hängt von der Masse der Zwerg- 
galaxie, dem Startradius r.; sowie von der ursprünglichen 
Exzentrizität ihres Orbits ab. 

Die Simulation wird in dieser Untersuchung genutzt, um 
die minimal nötige Masse rrimin der Zwerggalaxie zu be- 
stimmen, welche erforderlich ist, um sie vom Startradius 
Ti Über dynamische Reibung auf den Radius rt — 14 kpc 
innerhalb der simulierten Zeit von 13.7-10^ Jahren zu brin- 
gen. Die ursprüngliche Exzentrizität des Orbits muss 
berücksichtigt werden, da mit steigender Exzentrizität 
die Zwerggalaxie auf ihrem Orbit in Bereiche höherer Dich- 
te gelangt, welche die dynamische Reibung begünstigen. 
Dies ist auch aus Abbildung |5] ersichtlich: Je größer Ci, 
desto dichter kommt die Zwerggalaxie auf ihrem Orbit 
dem galaktischen Zentrum und daher den Bereichen ho- 
her Dichte. In Abbildung [TT] ist die nötige Mindestmasse 
fTT-rnin in Abhängigkeit des Startradius für verschiedene 
ei aufgetragen. 

Mit dem Startradius steigt die Mindestmasse rrimm- 
Bei Ci — ist zu erkennen, dass der maximale Startradius 
Vi im Bereich von 150 kpc liegt, da für größere Radien die 
Masse der Zwerggalaxie die des beobachteten Rings dunk- 
ler Materie überstiege. Bei kleinen Startradien und höher- 
en Exzentrizitäten reichen allerdings bereits Zwerggala- 
xien mit Massen um Faktor 5 kleiner als der Ringmasse, 
um bis zu rt innerhalb der simulierten Zeit zu gelangen. 
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Abbildung 11: Mindestmasse rrimin einer Zwerggalaxie, um 
vom Startradius r,; über dynamische Reibung innerhalb der si- 
mulierten Zeit rt zu erreichen. Der Zusammenhang von rrimin 
und Vi ist für verschiedene Exzentrizitäten dargestellt. 
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Tabelle 2: Parameter der Simulation des Einfalls mehrerer Ga- 
laxien. Masse m in Mq, Durchmesser d und Startradius in 
kpc, Inklination ii in Grad und numerische Exzentrizität e^. 
Die Durchmesser d werden so gewählt, dass die Zwerggalaxie 
um den Radius rt ~ lAkpc durch Gezeitenkräfte zerrissen wird, 
siehe Gleichung llSI Diese liegen im Bereich der für Zwerggalaxi- 
en dieser Größenordnung verwendeten und durch Beobachtung 
bestätigten Durchmesser, siehe zum Beispiel iBekkil (200^ ). 



3.4 Einfall mehrerer Galaxien 

Wie in Abschnitt 13.31 gezeigt, können auch Zwerggalaxi- 
en mit deutlich kleineren Massen als der Ringmasse bis 
zum Abstand von 14 kpc zum galaktischen Zentrum durch 
dynamische Reibung abgebremst werden. Dies ermögUcht 
die Vermutung, dass der Ring nicht aus einer großen, son- 
dern aus mehreren kleinen Zwerggalaxien entstanden sein 
könnte. An einem Beispiel soll nun konkret gezeigt wer- 
den, dass der Einfall mehrerer Zwerggalaxien ebenfalls zu 
einer ringförmigen Ansammlung kalter Dunkler Materie 
bei r = 14 kpc führen kann. 

Dabei wird der Einfall von drei Zwerggalaxien G„ mit 
Massen m zwischen 6 • lO^M© und 1.3 • IO^^Mq simuliert. 
Die Parameter der Galaxien sind in Tabelle [2] aufgelistet. 
Bei mehreren einfallenden Zwerggalaxienerhöhen jene, die 
den jeweiligen Abstand zum galaktischen Zentrum rt errei- 
chen, an welchem sie durch Gezeitenkräfte zerrissen wer- 
den, die relative Dichte der Milchstraße um diesen Radius 
(siehe Abschnitt 12.5p . Durch die gestiegene Dichte bei rt 
wird das Zerreißen einer weiteren einfallenden Zwerggala- 
xie an diesem Radius begünstigt. Dieser Prozess erschwert, 
dass Zwerggalaxien geringer Dichte weiter als bis zum Ra- 
dius Tt in die Milchstraße hinein geraten und trägt dabei 
sich selbst mit jeder zerrissenen Zwerggalaxie verstärkend 
zu einer Ringbildung um den Radius rt bei. Damit wird 
deutlich, dass der Einfall mehrerer Zwerggalaxien nicht 



zwingend mehrere Ringe mit sich bringt sondern, wie be- 
obachtet, das Zerreißen der Zwerggalaxien an einem be- 
stimmten Radius und damit die Bildung eines einzelnen 
Rings begünstigt wird. 

3.5 Rotationskurve 



In Abbildung [12] sind zwei Rotationskurven dargestellt. 
Die linke Rotationskurve ergibt sich aus der Simulation des 
Einfalls nur einer den Ring kalter Dunkler Materie verur- 
sachenden Zwerggalaxie mit der Masse m = 2.8 • IO^^Mq, 
dem Durchmesser d = 6.48kpc und einem Orbit mit 6^=0 
und ii = 0. Die rechte Rotationskurve stellt die Rotations- 
kurve der Simulation des Einfalls dreier Zwerggalaxien dar, 
deren Eigenschaften i n Tabelle [H aufgelistet si nd. Als Ver- 
gleichswerte sind von iHonma fc So"fu3 (Il997al) zusammen- 
gestellte Messwerte für w© = 220km/s und w© = 200km/s 
eingetragen, wobei Vq die Rotationsgeschwindigkeit der 
Sonne um das galaktische Zentrum bei r© = 8.3kpc be- 
schreibt. Die Fehlerbalken der Messung mit Vq = 200km/s 
entsprechen im Wesentlichen denen der Messwerte mit 
Vq = 220km/s und sind der Ubersicht halber nicht mit 
eingezeichnet. 

Befindet sich ein Körper zwischen dem CDM-Ring und 
dem galaktischen Zentrum, wirkt auf diesen eine Kraft in 
Richtung des galaktischen Zentrums, welche sich aus dem 
Gravitationspotential der sphärisch symmerischen Dichte- 
verteilung berechnen lässt, siehe Gleichung |ül In entgegen- 
gesetzte Richtung wirkt auf den Körper eine Kraft, welche 
sich aus dem Gravitationspotential des CDM- Rings ergibt, 
siehe Gleichung [5] Deshalb können Körper zwischen Ring 
und galaktischem Zentrum langsamer rotieren als sie es 
ohne den Ring könnten. Dieser Einfluss des Rings ist in 
den negativen Geschwindigkeiten in Abbildung [T^ ausge- 
drückt. 

Sowohl die Rotationskurve gemäß dem Modell des Ein- 
zeleinfalls einer Zwerggalaxie, als auch jene nach dem 
Modell des Einfalls mehrerer Z werggalaxien liegen inner - 
halb der Fehlerbalken der von iHonma fc So"fu3 (Il997ah 
zusammengestellten Messwerte der Rotationsgeschwindig- 
keit. Beide Modelle können den leichten Einbruch der Ro- 
tationsgeschwindigkeit bei r wlOkpc reproduzieren, stei- 
gen jedoch darauf etwas flacher an als die gemessenen 
Rotationsgeschwindigkeiten. Auch für große r liegen die 
Simulationsergebnisse noch innerhalb der Fehlerbalken. 
Letztere sind jedoch deutlich größer als im Bereich klei- 
ner Abstände zum galaktischen Zentrum. 



4 Diskussion 

An dieser Stelle werden die in Abschnitt [3] vorgestellten 
Ergebnisse der Simulation auf Bedeutung und Aussage- 
kraft hin überprüft. Dazu werden zuerst die Aussagekraft 
der Ergebnisse einschränkende Fehlerquellen bewertet. 

4.1 Fehlerquellen 

Die Ergebnisse entstammen einem semianalytischen Mo- 
dell des Einfalls einer und mehrerer Zwerggalaxien in die 
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Abbildung 12: Rotationskurve der Milchstraße mit Ring dunkler Materie als Ergebnis des Einfalls von einer (links) und drei 
(rechts) Zwerggalaxien. Die Beiträge der einzelnen Komponenten der Milchstraße zur Ge samtgeschwindigkeit sind dargestellt. 
Als Vergleichswerte (Messreihen mit vq — 220km/s und vq = 200km/s) sind von iHonma fc Sofue ( l_ 997a ) zusammengestellte 
Messwerte eingetragen. Beide Rotationskurven liegen im Wesentlichen innerhalb der Fehlerbalken dieser Messwerte. 



Milchstraße. Bezeichnend für ein semianalytisches Modell 
ist, dass sich die Ergebnisse lediglich aus den explizit mo- 
dellierten Prozessen und derer Zusammenwirkung erge- 
ben. Wichtig für aussagekräftige Ergebnisse ist also eine 
umfassende Vorstellung aller für das zu simulierende Sys- 
tem entscheidenden Prozesse. Modelliert wurde die gra- 
vitative Wechselwirkung zwischen der Zwerggalaxie und 
der Milchstraße mit durch Beobacht ungsdaten gest i itzten 



Dichteverteilungen der Milchstraße (jNavarro et al.l . 11996 



Dehnen et äD . ll998l) 



Das verwendete Gravitationspotential (siehe Abschnitt 
stellt zwar eine Näherungslösung dar, führt allerdings 
zu Rotationskurven innerhalb der Fehlerbalken von Mes- 



sungen der Rotationsgeschwindigkeit (jHonma fc Sofud . 
1997ah . Dies spricht für die Zulässigkeit der Näherungen. 

Das Modell zum Zerreißen der Zwerggalaxie durch Ge- 
zeitenkräfte berücksichtigt weder den Massenverlust der 
Zwerggalaxie auf ihrem Orbit, noch beschreibt es den Pro- 
zess des Zerreißens. Das Berechnen des Abstands rt zum 
galaktischen Zentrum, an welchem die Zwerggalaxie durch 
Gezeitenkräfte zerrissen wird, ist jedoch am Zer reißen der 
großen Magellanschen Wolke getestet. Der von 'Tremaine' 
( 1976f ) publizierte Radius rt ist knapp Faktor 2 größer als 
der selbst berechnete, allerdings berücksichtig t das hier 
verwe ndete Modell im Gegensatz zu dem von 'Tremaind 
( 1976h verwendeten die Auswirkungen des Gravitations- 
potentials der Disk. 

Die Reichweite der dynamische Reibung wird über den 
Coulomb-Logarithmus InA festgelegt. Dieser wurde mit 
Hilfe einer Computersimulation abgeschätzt, siehe Ab- 
schnitt 12.31 Je nach Geometrie des Halos und Position 
der Zwerggalaxie im Halo ändert sich mit dem für die 
dynamische Reibung wirksamen Bereichs auch InA, in 
diesem Modell wird jedoch vereinfachend ein konstanter 
Coulomb-Logarithmus verwendet. Dieser konstante InA 



liegt im Bereich der von Just &: Penarru bia ('2005') publi- 
zierten Coulomb-Logarithmen für das Milchstraßenhalo. 

Die durch Kollision abgebremste baryonische Materie 
der Zwerggalaxie befindet sich gemeinsam mit der CDM 
im Potential der Zwerggalaxie und führt daher durch gra- 
vitative Wechselwirkung auch zu einer Abbremsung der 
CDM. Allerdings ist die Masse der baryonischen Materie 
klein im Vergleich zu der Masse der CDM, weshalb die- 
se Form der Abbremsung der CDM nicht berücksichtigt 
wurde. 

Das Verfahren zur numerischen Zeitintegration des Or- 
bits wurde auf numerische Konvergenz getestet. Dazu 
wurden Orbits mit einer zwei Größenordnungen kleineren 
Schrittweite At als der sonst verwendeten simuliert. Es 
ließen sich numerisch keine Unterschiede feststellen. 

4.2 Bedeutung der Ergebnisse 

Die Simulationen können für physikalische Eigenschaften 
der Zwerggalaxie und ihres Orbits Ausschlussgrenzen set- 
zen. Die maximale anfängliche Exzentrizität des Orbits 
der Zwerggalaxie kann auf einen Wert von »0.2 begrenzt 
werden, siehe Abschnitt 13.11 

Die notwendige Masse der Zwerggalaxie, um vom Star- 
tradius ri bis zum Ringradius bei i? = 14 kpc über dyna- 
mische Reibung zu gelangen, ist für < 150 kpc deutlich 
geringer als die Ringmasse, siehe Abbildung [TT] Dies führt 
zur Annahme, dass der CDM-Ring möglicherweise nicht 
aus einem Einzeleinfall entstanden ist, sondern ebenso gut 
aus dem Einfall mehrerer Zwerggalaxien herrühren könnte. 
Sowohl das Modell des Einzeleinfalls als auch ein Modell 
des Einfalls dreier Zwerggalaxien kann die Rotationskurve 
der Milchstraße gut reproduzieren, siehe Abbildung [T^l 

Ergebnisse wie die eben vorgestellten wurden bereits in 
anderen Arbeiten über andere Methoden und andere Mo- 



10 



Raphael Errani: Modell zur Entstehung eines Rings aus Dunkler Materie in der Milchstraßenebene 



delle versucht zu erhalten. Es gibt n- Teilchen Simulationen 
und semianalytische Modelle, welche die Entstehung ei- 
nes CDM-Rings durch einen Einzeleinfa ll einer Zwerggala- 
xie mo delli eren, vergleiche zum Beis piel iKazantzidis et al. 



20081) und 'Penarrubia et al. (120051). Auch diese Modelle 



gehen von Orbits kleiner Inklination und kleiner numeri- 
scher Exzentrizität aus. Allerdings wird in diesen Modellen 
nicht auf die mögliche Entstehung des Rings aus mehreren 
einfallenden Zwerggalaxien eingegangen. 

Das in dieser Arbeit vorgestellte semianalytische Mo- 
dell ist also in der Lage, Aussagen über die Entstehung des 
CDM-Rings bei i? = 14 kpc in der Milchstraße zu machen, 
die sich mit den Ergebnissen von n-Teilchen-Simualtionen 
und anderer semianalytischer Modelle decken und die- 
se, trotz der Einfachheit des hier vorgestellten Modells, 
ergänzen. 

4.3 Bezug zu anderen Galaxien 

Auch in Rotationskurven anderer Galaxien finden sich, 
wie bei der Milchstraße, Einbrüche der Rotationsge- 
schwindigkeit bei bestimmten Abständen zum galakti- 
schen Zentrum. So zeigt NGC 660 (Hubble-Typ Sc) bei 
einem Abstand von etwa lOkpc zum galaktischen Zentrum 
einen Einbruch der Rotationsgeschwindigkeit, welcher ver- 
bunden mit einem Ring aus Sternen g roßer Inklination 
bezügl ich der galaktischen Ebene ist (jHonma fc Sofud . 

Auch eine solche Struktur lässt sich mit dem hier 
vorgestellten Modell erklären, setzt jedoch einen Orbit der 
eingefallenen Zwerggalaxie hoher Inklination voraus. Auch 
NGC 1808 (Sbc) und NGC 4258 (Sb) zeigen niedrigere Ro- 
tationsgeschwindigkeit um einen Abstand von lOkpc zum 
galaktischen Zentrum, allerdings wird kein Sternring bei 
diesem Radius beobachtet. Auch hier können eingefallene 
Zwerggalaxien für die Einbrüche der Rotationsgeschwin- 
digkeit verantwortlich sein, bei diesen Galaxien allerdings 
mit Orbits geringer Inklination, so dass die Sternströme 
innerhalb der galaktischen Ebene liegen. Das Finden ähn- 
licher Strukturen bei verschiedenen Galaxien deutet auf ei- 
ne ähnliche Entstehungsgeschichte und bekräftigt das hier 
vorgestellte Modell nach dem Prinzip der hierarchischen 
Galaxienentstehung. 

5 Zusammenfassung 

Um die Bildung eines Rings aus kalter Dunkler Mate- 
rie (CDM) innerhalb der Ebene der Milchstraße bei ei- 
nem Abstand von i? = 14 kpc zum galaktischen Zentrum 
zu erklären, wurde ein semianalytisches Computermodell 
entwickelt. Dieses Modell beschreibt die Entstehung des 
CDM-Rings gemäß dem Prinzip der hierarchischen Gala- 
xienentstehung durch den Einfall von Zwerggalaxien in die 
Milchstraße. 

Eine Zwerggalaxie befindet sich zu Beginn der Simu- 
lation bei einem Abstand zum galaktischen Zentrum 
der Milchstraße und besitzt ein Orbit mit der Anfangsge- 
schwindigkeit Vi, der Inklination ii und der numerischen 
Exzentrizität e^. Durch dynamische Reibung verliert die 
Zwerggalaxie an Drehimpuls und Bewegungsenergie und 



bewegt sich im Gravitationspotential der Milchstraße in 
Richtung des galaktischen Zentrums. Erreicht die Zwerg- 
galaxie den Abstand rt zum galaktischen Zentrum, über- 
steigen die durch das Gravitationspotential der Milchstra- 
ße wirksamen Gezeitenkräfte die Kräfte, welche die Zwerg- 
galaxie zusammen halten. Die Zwerggalaxie wird daher 
am Radius rt zu einem Sternstrom (bestehend sowohl aus 
Sternen als auch aus CDM und Wasserstoffgas) zerris- 
sen, der aufgrund der Ringform und der auf ein größe- 
res Gebiet verteilten Materie unanfällig für dynamische 
Reibung des Kalos ist. Daher behält der Sternstrom sei- 
nen Abstand zum galaktischen Zentrum in etwa bei. Ist 
die Inklination des Orbits der eingefallenen Zwerggalaxie 
klein, so kann wie beschrieben die Bildung eines Rings aus 
CDM innerhalb der Ebene der Milchstraße erklärt wer- 
den. Bei größeren Inklinationen können sich Sternströme 
außerhalb der galaktischen Ebene bilden. Das vorgestellte 
Modell ermöglicht außerdem. Aussagen über die physika- 
lischen Eigenschaften der Zwerggalaxie sowie ihres Orbits 
zu erhalten. Um in einer Zeit von maximal 13.7 • 10^ Jah- 
ren einen Ring aus CDM innerhalb der galaktischen Schei- 
be und mit geringer Exzentrizität zu erhalten, kann, so 
zeigen die Simulationen, die anfängliche Exzentrizität 
des Orbits der eingefallenen Zwerggalaxie nicht größer als 
~ 0.2 gewesen sein. Bei einer Masse des CDM-Rings von 
9 • 10^" Sonnenmassen ist die Mindestmasse, um von klei- 
nen ri < 100 kpc bis zum Radius von i? = 14 kpc zu gelan- 
gen, deutlich geringer als die Ringmasse. Dies ermöglicht 
die Annahme, dass der CDM-Ring nicht aus einer großen 
sondern aus mehreren kleinen Zwerggalaxien entstanden 
ist. Die Computersimulationen zeigen, dass sowohl ein Ein- 
zeleinfall als auch der Einfall mehrerer Zwerggalaxien die 
beobachtete Rotationskurve der Milchstraße gut reprodu- 
zieren kann. 

6 Fazit 

Das in dieser Arbeit vorgestellte semianalytische Compu- 
termodell ermöglicht Aussagen über die physikalischen Ei- 
genschaften von Zwerggalaxien, die einen Ring aus kalter 
Dunkler Materie in der Milchstraße nach ihrer Zerstörung 
durch Gezeitenkräfte gebildet haben. Die Ergebnisse des 
vorgestellten Modells decken sich mit denen anderer Publi- 
kationen und ergänzen diese um die Möglichkeit der Ring- 
bildung durch den Einfall mehrerer kleiner Zwerggalaxien. 
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